Unrestricted Hartree-Fock for Quantum Dots 
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We present detailed results of Unrestricted Hartree-Fock (UHF) calculations for up to eight elec- 
trons in a parabolic quantum dot. The UHF energies are shown to provide rather accurate estimates 
of the ground-state energy in the entire range of parameters from high densities with shell model 
characteristics to low densities with Wigner molecule features. To elucidate the significance of 
breaking the rotational symmetry, we compare Restricted Hartree-Fock (RHF) and UHF. While 
UHF symmetry breaking admits lower ground-state energies, misconceptions in the interpretation 
of UHF densities are pointed out. An analysis of the orbital energies shows that for very strong 
interaction the UHF Hamiltonian is equivalent to a tight-binding Hamiltonian. This explains why 
the UHF energies become nearly spin independent in this regime while the RHF energies do not. 
The UHF densities display an even-odd effect which is related to the angular momentum of the 
Wigner molecule. In a weak transversal magnetic field this even-odd effect disappears. 

PACS numbers: 73.21.La,31.15.Ne,71.10.Hf 



I. INTRODUCTION 

In the present work we discuss properties, predictions, 
and limitations of Hartree-Fock (HF) calculations for 
quantum dots. This method has a long tradition in 
atomic and nuclear physics, its application to quantum 
dots is therefore natural and has been discussed in var- 



ious recent papers 
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As we will demonstrate, 



some of the conclusions drawn on the basis of HF cal- 
culations are not based on firm grounds. This is in par- 
ticular the case, when the HF wave functions are used 
to describe charge distributions in a quantum dot. On 
the other hand, Unrestricted Hartree-Fock (UHF) will 
be shown to give rather reliable estimates for the ground 
state energies. 

While quantum dots may be considered as tunable ar- 
tificial atoms, the electron density can be much smaller 
than in real atoms and correlations play a more promi- 
nent role. 9 This is why for quantum dots the HF method 
has to be regarded with care. In this work we focus 
on the crossover from weak to strong Coulomb interac- 
tion, i.e. from higher to lower electronic densities. This 
is equivalent to weakening the external confinement po- 
tential for a given host material of the quantum dot. 

The physics of this crossover can be sketched as fol- 
lows: In the case of weak interaction (high density) a 
one-particle picture is valid: Electrons arc filled into the 
energy shells of the two dimensional isotropic harmonic 
oscillator. Here, the appropriate method is Restricted 
Hartree-Fock (RHF)ji^ where every orbital belongs to 
an energetic shell and has good orbital momentum. This 
shell filling with Hund's rule has been probed experimen- 
tally in small dots— In the case of strong interaction 
(low density) one can no longer stay within this simple 
one-particle picture: Wigner— has shown that for strong 
correlation the ground state of the 2D electron gas is 
described by localized electrons, representing a classical 
hexagonal crystal. Accordingly, in this limit the elec- 
trons in the dot form a small crystal, a so-called Wigner 



molecule, and the picture of energetic shells is no longer 
meaningful. One has to improve the HF approximation 
by passing over to UHF which means that the space of 
the HF trial wave functions is extended. The UHF Slater 
determinant lowers the energy by breaking the symmetry 
of the problem, i.e. spatial and spin rotational invariance. 
This complicates the interpretation of the UHF solution. 

For very strong interaction UHF is also expected to 
give reasonable results because a one-particle picture of 
localized orbitals^ should model the Wigner molecule 
quite well. In fact, the UHF energies become nearly 
spin independent, while this is not the case with RHF. 
We show that the UHF Hamiltonian for strong interac- 
tion has the same spectrum as a tight-binding Hamilto- 
nian of a particle hopping between the sites of a Wigner 
molecule. The hopping matrix elements and on-site 
energies can be extracted from the UHF orbital ener- 
gies. The localization-delocalization transition has al- 
ready been probed experimentally in larger quantum 
dots^ so Wigner molecule spectroscopy is within reach 
of current technology. 

An incomplete account of our results has been pre- 
sented in an earlier short communication^ Here, we dis- 
cuss in detail the two-electron problem and present an 
elaborate analysis of the limit of strong interaction. In 
Sect. |H] we shortly recall the model and method. In 
Sect. ITTT1 we obtain explicit results for quantum-dot He- 
lium that already show many features of HF solutions for 
higher electron numbers presented in Sect. II VI In Sect. El 
we also discuss the effect of a magnetic field. 



II. HAMILTONIAN AND HARTREE-FOCK 
APPROXIMATION 

In this work we follow the notation and method pre- 
sented in our earlier article^ for zero magnetic field. The 
Hamiltonian of an isotropic parab olic qua ntu m dot with 
magnetic field reads (see e.g. Refs. 
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where the positions (momenta) of the electrons are de- 
noted by rj (pj). The effective mass is m* , and the 
dielectric constant is k. The vector potential of a ho- 
mogeneous magnetic field B orthogonal to the plane 
of the quantum dot in symmetric gauge reads A(r) — 
V) x, 0), and the corresponding cyclotron frequency 
is u c — eB/m* . 

Now we can introduce oscillator units, and describe 
the system di mensionless: energies in units of fru> c fi = 
hy/uj 2 + w 2 /4 and lengths in units of l = ^h/m*uj c s- 
Then the Hamiltonian takes the form 

where we have introduced the dimensionless coupling 
constant 

\ = Io/cib = s / nlofiuj (3) 

with the effective Bohr radius a B . For example A = 2 
corresponds to fuwSmeV for a GaAs quantum dot. The 
Hamiltonian J5J is formally the same as without mag- 
netic field, apart from an additional term proportional 
to the total angular momentum which scales with the di- 
mensionless parameter— f2 c := w c /co e g. The major part 
of our calculations presented below is for zero magnetic 
field. 

Regarding the HF approximation, 25 let us recall the 
expansion of the HF orbitals in terms of the angular mo- 
mentum eigenfunctions of the two-dimensional harmonic 
oscillator 6 

(r|i> = Vi (r) = Yl <M(r\nMai). (4) 

n=0,oo 
M= — oo,oo 

Here, M is the angular and n the radial quantum number 
of the Fock-Darwin basis. Each orbital has its own fixed 
spin <Ji — ±1/2, this means there is no double occupancy 
of orbitals with spin up and down, but there are different 
orbitals for different spins. Thus only the z-component of 
the total spin is fixed, S 1 * * = J2i &i = S z . Furthermore, 
the orbitals Q are in general no longer eigenfunctions 
of the one-particle angular momentum (UHF). Therefore 
the HF Slater determinant is not an eigenstate of the 
total angular momentum L^ ot , it breaks the symmetry 
of the original Hamiltonian. 27 Another possibility is to 
give each orbital i a fixed angular momentum Mi. With 
this restriction one obtains RHFi^ which preserves the 
total angular momentum but yields higher ground-state 
energies. Still another possibility is to build a Slater de- 
terminant of spatially localized orbitals for the strongly 
interacting case£ or of multicenter localized orbitals in 



high magnetic field 8 and vary these orbitals to minimize 
the HF energy. Our orbitals are self-consistent and are 
best adapted to study the crossover from weak to strong 
correlation. 

In principle the orientation of the deformed symmetry- 
breaking HF solution is arbitrary. This is due to the ro- 
tational invariance of the original Hamiltonian and can 
be called orientational degeneracy. The actual UHF so- 
lution found has a special orientation and it depends on 
the initial guess for the density matrix. Often but not 
always the symmetry breaking is manifested in the HF 
single-particle density n HF (r) = X)i=i l^ 7 *)! 2 - For a 
quantum dot in zero magnetic field, the Hamiltonian is 
invariant under time reversal. Thus we can choose real 
expansion coefficients u l nM in However, then the 

HF one-particle density is always symmetric to one axis. 
Any arbitrary orientation can be obtained by applying 
exp(iaL* ot ) to the Slater determinant. 

III. UNRESTRICTED HARTREE-FOCK FOR 
QUANTUM-DOT HELIUM 

In this section we present UHF energies and densities 
for the two-electron quantum dot (quantum-dot Helium) 
at zero magnetic field for increasing interaction strength 
A. This illustrates the basic concepts and properties of 
the HF approximation, and reveals features that are also 
important for higher electron numbers. We compare with 
exact results obtained by diagonalization of the relative 
motion. We also compare with the RHF method, in order 
to illustrate the differences to UHF. 

The UHF two-electron problem has been treated pre- 
viously by Yannouleas and Landman. 5 However, we find 
some deviations from their results. An extensive dis- 
cussion of the RHF solution for quantum-dot Helium at 
A ~ 2 can be found in Rcf. [H Finally, we want to men- 
tion that the two-electron problem has also an analytic 
solution in terms of a power series— 

A. Two-electron Slater determinant 

The Slater determinant for two electrons with S z = 

is 

= 7! M r iV2(r 2 )x+X- " Vi{r2)m{ri)X 2 +X-} ■ 

(5) 

Here we have displayed the orbital and spin parts of the 
wave function explicitly, x± 1S f ne spin of the z-th elec- 
tron. The state \E ,HF is generally not an eigenstate of the 
total spin S 2 ot . In order to obtain a singlet one has to 
set tpi = (p2, and thus 

* HF = ^l(nVl(r2)Xsinglet- (6) 

This restriction is also called closed-shell HF (CSHF), 
because if every orbital is filled with spin up and down, 
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FIG. 1: Comparison of different HF energies for quantum-dot 
Helium with the exact ground state energy vs. the coupling 
constant A. 



open shells are impossible. One sees from iJSJ that the 
Slater determinant violates the symmetry of the problem. 
For two electrons the spin symmetry is easily restored, 
namely by a superposition of two Slater determinants 
with spin up/down and down/up. For the polarized case 
S z = 1, the total spin is conserved, and the HF wave 
function is a product of a symmetric spin function and 
an antisymmetric orbital function. 




■0-0 S =1 exact 
S 2 =0 UHF 

B£]S z =lUHF 

q..q S=ORHF 

»-oS_=lRHF 



FIG. 2: Absolute energy differences with respect to the exact 



S = ground state AE = E s - Sgs" 
UHF energies are nearly the same. 



Above A « 4 the two 



within about 0.3%, the S z = 1 state is somewhat lower 
than the S z = state. The exact energies merge more 
slowly: for A = 20 the energy difference between singlet 
and triplet is still about 1%. Note that the RHF energies 
fail to become spin independent for large A, as can be 
seen from Fig. [3 Of course, one expects spin indepen- 
dent energies in the classical limit of localized electrons 
without overlap. 



C. UHF one- particle densities 



B. Different HF approximations 

We now compare the energies of different HF approx- 
imations with the results of an exact diagonalization. 6 
First we consider the case S z = 0. The most general 
ansatz for the HF orbitals is the UHF state J2J, a spin 
dependent expansion with arbitrary angular momentum. 
Less general is the RHF ansatz, where angular momen- 
tum is preserved. And still less general is CSHF JBJ, 
when we force the two electrons to occupy two identi- 
cal (rotationally symmetric) orbitals. In Fig. ^ on e can 
clearly see the importance of breaking the symmetry to 
obtain lower HF energies. Up to Aw 1 all three methods 
give nearly the same result. Up to A ~ 3 the closed- 
shell energy is equal to the RHF energy. In other words: 
From this point on the two RHF orbitals are no longer 
identical. As expected the UHF energy is lowest. 

In Fig. [3 we show the differences of the RHF and UHF 
energies from the energy of the exact ground state which 
is the singlet. For S z = 1 one needs two different or- 
bitals, there is no CSHF. The UHF method gives lower 
energies than RHF, but the gain in energy is not as big as 
in the unpolarized case. Interestingly, the UHF energies 
become spin independent with increasing A: they agree 



Now we want to have a closer look at the one-particle 
density which is just the sum of the densities of the two 
orbitals, n HF (r) = \(fii(r)\ 2 + |^ 2 (t")| 2 . 

In Figs. 01 and 0] we show this density for different 
values of the coupling parameter A. Already for a rel- 
atively small A we detect two azimuthal maxima. The 
density is strongly anisotropic which is due to the sym- 
metry breaking. In the case of S z = 1 the two maxima 
are more distinct as a consequence of the Pauli prin- 
ciple: spin-polarized electrons are more strongly corre- 
lated. However, the direct interpretation of the two dips 
as localized electrons is questionable. With increasing A 
the azimuthal modulation first decreases, but for A > 8 ( 
A > 10 for S z — 1) it increases again. For very high A the 
densities become almost spin independent. A closer view 
reveals that the azimuthal maxima are more distinct for 
the case S z = 0. This arises from the exchange term in 
the HF energy: it lowers the energy for strong interaction 
and overlapping spin-polarized orbitals. 

While the azimuthal modulation is an artifact of the 
HF method, the densities display correctly a minimum 
in the center which gets deeper with stronger interaction. 
Also, the maxima are in very good agreement with the 
classical positions r a — ^/A/4 (see Appendix ^J. 
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FIG. 3: Shadowed contour plots of the UHF one-particle den- 
sities n HF for TV = 2, S z = 0. One contour corresponds to 
1/10 of the maximal density, (a) A = 2, (b) A = 6, (c) A = 8, 
(d) A = 20. 
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FIG. 5: Pairs of real UHF orbitals for N = 2, S z = 1. (a), 
(b) at A = 2, (c), (d) at A = 10. 



D. UHF orbitals and orbital energies 
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FIG. 4: UHF one-particle densities for N = 2, S z = 1. (a) 
A = 2, (b) A = 6, (c) A = 8, (d) A = 20. 



In order to understand the form of the UHF densities 
it is useful to have a closer look at the UHF orbitals. 
For S z = we find two orbitals that are exactly complex 
conjugate, ip\ = ip^. This can be seen by studying the 
expansion coefficients u n M in Eq. Q and means that the 
Slater determinant is symmetric under time reversal. 

For S z = 1 the two orbitals depicted in Fig. [S] are 
always different and can be chosen real. For A = 2 one 
can still interpret the orbitals in the energy shell picture 
of RHF: the first orbital is (approximately) round, S-likc, 
and the second one is dumbbell formed, P-likei2£ 

For very high A > 14 there is a simple relation between 
the orbitals for the two spin polarizations: for S z = 1 we 
may choose both orbitals real and then we find 

^c^i^r 1 ). (7) 

In this fashion, we see that ffj^ are complex conjugate 
and approximately orthonormal. 

To shed more light on this behavior we consider also 
the orbital energies. We start with the HF Hamiltonian 
in the HF basis for S z = 1 

( Ei \ / hn+w 12 i2 \ , . 

\0 £2 ) \ h 22 +w 12 i2 ) ' [ ' 

Here, we use the notation hij = (i\h\j) and ify/y = 
(ij\w\kl) for matrix elements in the HF basis (see Ref.0). 
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FIG. 6: UHF Conditional probability density n HF (x\y) for 
TV = 2, S z = 1. In the upper row y = (2, 0) (x) (a) A = 2, (b) 
A = 6. Lower row y = (0, 2) (x) (c) A = 2, (d) A = 10. 



When wc apply the unitary transform l|7|l 



Hi 



1/1 i\ e x Q\ 1 1\ _ U -t 
2 v. 1 -* J \ e 2 J \ -i % ) \-t U 

we obtain a two-state Hamiltonian H2, with on-site en- 
ergy U = (/in + h?,2 + 2 , u;i2i2)/2 and tunnel splitting 
t = Q122 — hu)/2. Thereby we have mapped the HF 
Hamiltonian on a lattice problem.lt is intuitive that for 
strong interaction the two electrons localize, and thus 
a tight-binding approach should become physically cor- 
rect. This is also the case for larger electron number as 
discussed below. 



E. UHF two-particle densities 

Next we examine the conditional probability density 
(CPD) for finding one electron at x, under the condition 
that another electron is at y. For quantum-dot Helium 
and S, = the CPD reads 



n Hb (x\y) = 



\<pi(*)\ 2 \<P2(v)\ 2 + \<pi(v)\ 2 \Mx)\* 



]HF 



(v) 



(10) 



Now, 



since we found complex conjugate orbitals, <pi 

HF,' |„,\ _ „HF( 



(P2, we have n (x\y) = n (x), i.e. the conditional 
probability density is independent of the condition. This 
is not really astonishing, because within the HF method 
two electrons are only correlated by the exchange term, 
which vanishes here^ 



For S z = 1 the orbitals are different from each other 
and the CPD is given by 

n HF (x\y) = {\ Vl {x)\ 2 \ V2 (y)\ 2 + \^(y)\ 2 \M^)\ 2 

~2ReK( a; )^2^)^i(y)^(y)]}/« HF (y) • (n) 

In Fig. El we show contour plots of UHF CPDs for dif- 
ferent coupling constants and given positions y. In the 
upper row, for y = (2,0), we find for small A = 2 a 
suggestive result: the density has a single maximum at 
a distinct distance from the fixed coordinate y. With 
increasing A, however, we obtain two maxima, which de- 
velop more and more and are not at all located at the 
classical position. 

The situation is likewise irritating when one chooses 
y = (0,2) as fixed coordinate (lower row). While the 
exact CPD is rotationally symmetric when both x and 
y are rotated, the UHF CPD does not respect this sym- 
metry. The reason for this lies in the symmetry breaking 
which cannot completely account for correlations. The 
UHF Slater determinant is deformed and derived quanti- 
ties do not necessarily have a direct physical meaning, - 
except for the UHF energy which is a true upper bound 
for the exact energy. 



IV. UNRESTRICTED HARTREE-FOCK FOR 
HIGHER ELECTRON NUMBERS 

In this section we show further results of UHF cal- 
culations, namely energies and densities for up to eight 
electrons (B — 0). Many effects are similar to what we 
have already seen for two electrons, for example the er- 
rors of the UHF energies and their spin dependence. An 
interesting phenomenon shown by the UHF densities is 
the even-odd effect discussed below. 



A. UHF energies 

For N > 2 we compare the UHF energies with results 
of a Quantum Monte Carlo (QMC) simulation by Egger 
et al^i These results were obtained for a very low tem- 
perature T = O.lhuj/kB- The QMC energies are always 
below the HF energies and can therefore be considered 
as effective zero temperature reference points. 

For N — 3 QMC, a semiclassical analysis^ as well 
as an exact diagonalization study^i predict a transition 
from the S = 1/2 ground state in the weakly interacting 
case to a S = 3/2 ground state for A > 4. Within UHF 
this transition occurs already near A = 2. In Fig. one 
can see that the relative error for S z = 3/2 is small, 
less than 3% . In the non-polarized case the error is 
higher, about 7% for A > 2. With increasing N and A 
the relative error becomes smaller because the absolute 
energies are higher. 

In Fig.[S]we show the absolute energy differences from 
the QMC ground state for eight electrons. For intermedi- 
ate values of A the UHF energies become already nearly 
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FIG. 7: Relative error of the UHF energy (E$ F 
E® MC )/Ef MC for N = 3 vs. coupling constant A. 
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FIG. 9: Even-odd effect of the UHF one-particle densities 
n HF for A = 6, different electron numbers N and polarized 
spin S z = N/2, (a) AT = 3, (b) N = 4, (c) AT = 5, (d) N = 8. 



UHF renders the properties of the spin-polarized solu- 
tion for larger A. This can also be seen in the UHF den- 
sities, which become spin independent for larger interac- 
tion (see below). Finally, we briefly mention the RHF 
results: there for large A the HF energies do not become 
spin independent, but the energies for lower spins are con- 
siderably higher. For large A RHF gives a poor estimate 
of the ground state energy. 



B. HF densities: Even-odd effect 



FIG. 8: Absolute energy differences from the QMC ground 
state, AE = Es — Eq S for eight electrons and various spins 
vs. coupling constant A. 

spin independent, whereas the QMC energies approach 
this semiclassical behavior more slowly. For stronger in- 
teraction the HF ground state is always spin-polarized. 
Thus the UHF method can not resolve the correct spin 
ordering of the energies. 

For N — 8 the QMC method predicts a crossover of 
the total spin from S — 1 to S — 2 near A = 4. The 
UHF method finds a polarized ground state with 5 = 4 
for A > 4. There, however, the energy differences for 
different spins are already quite small. 

One can conclude that the UHF Slater determinant 
with fixed spin structure gives a rather poor description 
of the total many-electron wave function. Essentially, 



In this subsection we consider the UHF densities for 
higher electron numbers. We first show in Fig. [5] the 
densities for rather strong coupling constant A = 6, var- 
ious electron numbers N and S z — N/2. Above this in- 
teraction strength the UHF densities are essentially the 
same for all S z (except for N = 2, see above) and do not 
change qualitatively with increasing A. 

Surprisingly, only for some N does one obtain a 
molecule-like structure, i.e. an azimuthal modulation as 
seen for two electrons. For three and five electrons the 
density is apparently rotationally symmetric and also for 
eight electrons, where we have a pronounced maximum in 
the center. The expected molecule-like structure shows 
up only for N — 2 and 4. Thus, when we consider also 
N = 6, 7 (see below) we recognize that azimuthal max- 
ima occur only for an even number of electrons per spatial 
shell. In stating this we want to emphasize, that all the 
densities shown belong to symmetry breaking, deformed 
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Slater determinants. 

This even-odd effect is also surprising, because UHF 
calculations for quantum dots in a strong magnetic field^ 
found molecule-like densities for all electron numbers, 
and frequently a magnetic field leads to similar effects 
as a stronger interaction. We also have performed calcu- 
lations with a magnetic field that reproduce the densities 
of Ref. |j and show that the molecule-like structure dis- 
appears for odd N for vanishing fields 

A physical explanation of the even-odd effect combines 
the geometry of the classical system with the symmetry of 
quantum mechanics,— Consider the exact spin-polarized 
iV-electron wave function for the Wigner molecule 
case. Due to the strong Coulomb repulsion, the electrons 
move on an TV-fold equilateral polygon (for N < 6; for 
N = 6 one electron enters the center of the dot). A 
rotation by 2tt/N therefore corresponds to a cyclic per- 
mutation 

expj^Lr} *n = (-If' 1 *iv , (12) 

where we have used that a cyclic permutation of an 
even (odd) number of electrons is odd (even). From 
Eq. I|12|) the allowed total angular momenta of the 
Wigner molecule can be easily read off: for an odd num- 
ber of electrons the minimal angular momentum is zero, 
whereas it is nonzero and degenerate for an even elec- 
tron number, e.g. M tot = ±2 for N — 4. Hence, the 
UHF wave functions for N — 2,4,7 can be interpreted 
as standing waves, i.e. superpositions of opposite angu- 
lar momentum states. For odd numbers of electrons in 
a spatial shell there is no angular momentum degener- 
acy and therefore no standing wave and no modulation 
in the densities. With a similar argument Hirose and 
Wmgreeni^ explain the charge-density- waves which they 
found for odd number of electrons in the weakly interact- 
ing regime from density functional calculations. 

Equation l|12|l does not hold anymore when the spins 
are not polarized, because the total wave function is not 
a product of spin and orbital wave functions. However, 
within UHF we do not fix the exact spin but only sub- 
spaces with fixed S z . For S z < N/2 and strong interac- 
tion the UHF solution mainly renders the properties of 
the spin-polarized solution, since the energies and den- 
sities are essentially the same for A > 6. The even-odd 
effect is thus not a physical effect but an artifact of the 
UHF symmetry breaking. Therefore great caution must 
be taken when interpreting the UHF densities. In par- 
ticular, the exact onset of Wigner crystallization cannot 
be determined reliably from UHF calculations. 



C. Closer look at three electrons 

As we have just discussed, for three electrons with 
strong interaction we do not find the naively expected 
density with three maxima but a nearly round density. 



When we plot the density of Fig. H2a) with more con- 
tour lines (not shown) a tiny sixfold modulation of the 
density is discernible. This can be understood by going 
back to Eq. after M tot = the next allowed total 

angular momentum values are M tot = ±3, which give 
rise to a standing wave with six maxima. This becomes 
also clear from the densities of the single orbitals building 
the UHF single-particle density. In Fig. ^| we show the 
orbital densities for A = 4 and A = 6. We find a sixfold 
orbital, as well as two diametrically oriented threefold 
orbitals. One clearly recognizes how the sixfold modula- 
tion results from this. Note that the HF orbitals are not 
localized (for example at the angles of a triangle). 

At this point we want to address a related issue, the 
uniqueness of the HF orbitals. One can easily show with 
the help of the HF equations that HF orbitals with the 
same spin are no longer unique, if the corresponding one- 
particle energies Si are degenerate. In this case, any uni- 
tary transformation of degenerate orbitals also fulfills the 
HF equations. In Fig. ^3 the energies Si are degenerate 
for the two states (b),(c) and (e),(f). Therefore these 
two orbitals are no longer uniquely determined, - in ad- 
dition to the orientational degeneracy of the total Slater 
determinant which is physically obvious. 

Now, we want to have a closer look on the orbital ener- 
gies: it is natural to presume that their degeneracies are 
a signature of Wigner crystallization, i.e. the geometry of 
the Wigner molecule. For strong interaction one should 
be able to represent the system as a lattice problem on 
an equilateral triangle. The corresponding Hamiltonian 
for N = 3,S Z = 3/2 reads 



u 


-t 


-t 


-t 


u 


-t 


-t 


-t 


u 



(13) 



where U is the on-site energy and t is the tunneling ma- 
trix element between localized states. The eigenvalues of 
H3 are £\ = U — 2t and twice £2/3 = U + t which is in fact 
the degeneracy of the UHF orbital energies (Fig. I1U|) . 

On the other hand, for S z = 1/2 the tight-binding 
Hamiltonian involves tunneling only between the two 
spin up states and takes the form 



H's = \-t U . (14) 
\ u J 

The eigenvalues are £1/2 = U ± t (spin up) and £3 = U 
(spin down). With UHF for A = 6 we find E\ = 6.65, 
£2 = 7.10 and £3 = 6.87, which has to be compared with 
the orbital energies for the polarized state given in Fig.llUI 
and yields t sa 0.22. For larger A the agreement becomes 
better, e.g. for A = 12 we find e x = 10.140, £ 2 = 10.309 
and £ 3 = 10.224 for S z = 1/2, while £l = 10.06 and 
£ 2/3 = 10.313 for S z = 3/2, which gives t w 0.084 in 
both cases. 
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FIG. 10: UHF orbital densities \ipi\ 2 (i = 1,2,3) for N = 3 and S z = 3/2. Upper row A = 4, lower row A = 6. For the 
single-particle energies we obtain (a) s\ = 4.92 and (b), (c) £2 = £3 = 5.84; (d) E\ = 6.44 and (e), (f) £2 = £3 = 7.11. 



D. Lattice Hamiltonian and localized orbitals 

For large A the HF Hamiltonian has the same eigen- 
values as a lattice Hamiltonian. Thus, there must be 
one-to-one correspondence between these two. Remem- 
ber, however, that HF is a one-particle picture and thus 
the tight binding Hamiltonian describes one particle hop- 
ping on a grid. The HF Hamiltonian is diagonal in the 
HF basis Q, 



A' 



(15) 



Now, if the eigenvalues Si coincide with those of a lattice 
Hamiltonian, e.g. H3 in this means that we have 

to transform the UHF orbitals with the inverse of the 
orthogonal transformation which diagonalizes the lattice 
Hamiltonian to pass over to localized orbitals. The Slater 
determinant is not changed when we transform among 
occupied orbitals^ 



N 



(16) 



In this new basis the HF equations read 

N N 

Y,i(p\ h i«> + E(hh<f)K = ■ ( 17 ) 

q r 

Now, in the basis \p), we should have non vanishing 
(p\ h \q) only for nearest neighbors 31 and the contribution 



of the two-particle matrix element should essentially be 
given by the direct term, i.e. diagonal elements of the 
Coulomb interaction. Then 117fl reduces to 



JV 



N 



(18) 



which is now of the form of a lattice Hamiltonian. 

We now present strong numerical evidence for this 
connection between the UHF Hamiltonian and a lattice 
Hamiltonian for N — 4 and 5 which are the simplest cases 
of electrons on a ring. For N = 4, S z = 2 we have 



U -t -t 
-t U -t 

-t u -t 
-t -t u 



(19) 



with the eigenvalues E\ = U — 2t, e 2 /3 = U and £4 = 
U + 2t. The eigenvectors of H4 determine the transfor- 
mation l|l(j|) . Applying this transformation to the HF 
Hamiltonian, as wc did in JJjJ, we obtain for A = 8 
an Hamiltonian of the form 119(1 with U = 10.924 and 
t = 0.195. The next nearest neighbor hopping matrix 
element (hopping along the diagonal of the square) is 
t* = 2e 2 — £1 ~ £4 = 0.003, which is indeed very small. 

Likewise we can determine the lattice Hamiltonians for 
other electron numbers and spin configurations and we 
have collected results for t and U for stronger interaction 
up to A = 20. For N = 4, S z = 1 the lattice Hamiltonian 
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reads 



0.60 



U -t 

-t U -t 

-t U 

£7 



(20) 



with the eigenvalues £1 = U — \/2t, £2 = U and £3 = 
U + \/2i (spin up) and £4 = [/ (spin down), while for 
iV = 4, 5 2 = we have 



u -t 

-t f7 

u -t 

-t u 



(21) 



with e\/2 = U ±t (spin up), £ 3 / 4 = U ±t (spin down). 
Here, we have to assume that the four states are occupied 
with two pairs of nearest neighbor parallel spins in order 
to obtain agreement with the UHF orbital energies. The 
values of t we obtain in this way for the three spin states 
S z = 0, 1, 2 agree within 1% for A = 8 . 

For N — 5 we have a pentagon and again three dif- 
ferent spin states. For S z = 5/2 the lattice Hamiltonian 
with nearest neighbor hopping is 



( u -t -t 
-t u -t 
-t u -t 
-t u -t 
\-t -t u J 



(22) 



with the eigenvalues e% = U — 2t, £ 2 / 3 = U + t(l — Vo)/2 
and £4/5 =U + t(l + \/5)/2, while for S z — 3/2 we have 



H' 



( 



\ 



u -t 

-t u -t 

-t u -t 

-t u 

u J 



(23) 



with £ 1/2 = U- t(VE± l)/2, £ 3/4 = U + t(\/Z =F l)/2 
(spin up) and £5 = U (spin down). Finally for S z = 1/2 
we have 







V 





c/-too 

-! [/ 
u 
00c// 



(24) 



with the eigenvalues £1 



/ 3 = [/ ± t, £ 2 = f/ (spin up) 
and £4/5 = {/ (spin down). Note that here the values of 
the UHF orbital energies suggest a model with only two 
nearest neighbor parallel spins. For A = 6 the values of t 
for all three spin states coincide within 1%. 

Figure 1111 summarizes our findings about the tunnel 
matrix elements. Reference Il6l predicts t oc exp(— •JrT), 
where r s is the nearest neighbor distance of the electrons 
measured in units of the effective Bohr radius. Since 



0.10 



o N=2 
□ N=3 
N=4 
a N=5 




0.02 



FIG. 11: Log-linear plot of tunnel matrix element t vs. A 2//3 
for various electron numbers. For A > 8 the line of best fit is 
shown. 



classically r s oc A 4 / 3 (cf . Appendix we plot In t ver- 
sus A 2 / 3 . For A > 8 we find indeed a linear behavior. 
For lower A, the tunneling matrix element is not really 
defined, since the lattice model is not appropriate. The 
tunneling matrix element is largest for N — 2 because 
two electrons are always closest (see Appendix^ . Three 
electrons always have the smallest value of t because the 
corresponding equilateral triangle has a longer side than 
the square and the pentagon. For higher electron num- 
bers one electron enters the center of the dot, and the 
UHF spectra are more complicate but still show the typ- 
ical degeneracies. However, now the lattice Hamiltonian 
has various tunneling constants and on-site energies. 



E. Seven-electron Wigner molecule 

Seven classical electrons form a equilateral hexagon 
with one central electron, which is a fragment of a hexag- 
onal lattice. In Fig. ^| we show UHF densities for 
N = 7 starting with a small A. The UHF ground state is 
S z = 1/2 up to A < 3, then spin-polarized. In Fig. lbW a) 
for A = 1 we see a fourfold modulated density. How is 
that possible for seven electrons? The answer is that in 
this case the energy shell picture of the harmonic oscil- 
lator is still valid: six electrons are just a shell closure 
and the next electron is put in the new shell in an orbital 
with maximal angular momentum. This angular momen- 
tum is M — ±2 and from the su perp osition one obtains a 
fourfold standing wave (cf. Ref . Il4|) . Here, the energy is 
basically the same as in RHF, but the Slater determinant 
breaks the symmetry. 

With increasing interaction strength a Wigner 
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-4 -3 -2 -1 1 2 3 4 -4 -3 -2 -1 1 2 3 4 



FIG. 12: One-particle densities for the UHF ground state of 
N = 7 electrons, (a) A = 1, (b) A = 2, both S z = 1/2. (c) 
A = 4, (d) A = 10, both S z = 7/2. 



molecule is formed with one electron in the center and 
six in the surrounding ring [Fig. H2f bWcl], We want 
to emphasize that the UHF densities mirror the classi- 
cal shell filling. This can even be quantified: the posi- 
tions of the maxima (even in the 'round' densities) agree 
very well with the classical configurations in Appendix 1X1 
From the UHF density the nearest neighbor distance f s 
can be determined. For example from Fig. ll2f d') we find 
r s ~ 3.0, which is also the classical value. Here we have 
to take into account that we measure length in oscillator 
units. Frequently, one is interested in the density param- 
eter r s given in effective Bohr radii. 32 Then Fig. I12f d~) 
gives r s = r s lo/a^ — Xf s « 30. The r s values we obtain 
in this way agree also well with the results of Ref. IT7I 
There r s is determined from the first maximum of the 
two-particle correlation function. 



V. UNRESTRICTED HARTREE-FOCK WITH A 
MAGNETIC FIELD 

In this section we want to present some calculations 
with a magnetic field orthogonal to the plane of the quan- 
tum dot. This system has been discussed extensively in 
the literature, especially in connection with the quantum 
Hall effect. UHF calculations by Miiller and Koonin 4 
have shown a magnetic field induced Wigner crystalliza- 
tion. However, they only considered the limiting case 
of a strong magnetic field and therefore included in the 
basis for expanding the UHF orbitals only states from 




-4 -3 -2 -1 1 2 3 4 -4 -3 -2 -1 1 2 3 4 

FIG. 13: Evolution of the UHF one-particle densities for 
AT = 3, S z = 3/2 and A = 10 with increasing magnetic field 
strength u> c — uj c /uj. (a) Q c = 0, (b) Q c — 0.5, (c) u) c = 1.5, 
(d) Q c = 2.5. 



the lowest Landau level (Fock-Darwin levels with n = 0). 
The high field case has also been studied by Palacios et 
al£ and Ruan et alpi£*i2*2S To study smaller magnetic 
fields, our basis is better adjusted to the problem. It is 
intuitively clear, that electrons are further localized by 
the magnetic field. Indeed, for sufficiently strong fields, 
we do not find an even-odd effect for UHF densities but 
molecule-like densities for all electron numbers. 

Numerically, thanks to the similar form of the Hamil- 
tonian J2J) to the one without magnetic field, the general- 
ization of our UHF code is straightforward. However, the 
magnetic field breaks time reversal symmetry, left and 
right turning solutions are no longer energetically degen- 
erate. Therefore in the expansion of the UHF orbitals 
Q we have to use complex coefficients. 

We first consider three electrons and a large interaction 
parameter A = 10. This means that we have a shallow 
quantum dot where the Coulomb interaction dominates 
and the magnetic field is relatively weak. In Fig. 1131 we 
display the evolution of the UHF one-particle densities 
with increasing magnetic field strength uj c — uj c /^> at 
fixed A. This is not exactly the physical situation, cor- 
responding to a quantum dot exposed to an increasing 
magnetic field, since the coupling constant A becomes 
smaller with increasing field. Here we just want to show 
that a magnetic field does not have the same effect on 
the UHF density as a strong interaction. 

In Fig. 113r d) we see three distinct, localized electrons 
in the UHF density. The three single orbital densities 
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have nearly the same form. They are thus similar to 
the orbitals chosen in Ref. 0. With decreasing magnetic 
field strength the maxima in azimuthal direction vanish 
slowly, until we have again a nearly round density for 
uj c = as in Fig.^a). The density in Fig. ll3f a) has been 
obtained from an initial guess with threefold symmetry. 
Therefore we can be sure that we have not obtained a 
local minimum but the true HF ground state. 

As a second example we show the evolution of the UHF 
density of six electrons at intermediate coupling strength 
A = 3.2. Without magnetic field the density is round, 
Fig. I14f a). and with a weak magnetic field fivefold with 
a central electron, Figs. I14f b').fc % ). Remarkably, for in- 
termediate magnetic field Q c w 1 . . . 1.5, the UHF ground 
state has a perfectly round density, Fig. I14f d). and also 
a rotationally symmetric Slater determinant. This is 
the so-called maximum-density-droplet of MacDonald et 
al.j2£ where the electrons occupy the lowest orbitals with 
increasing angular momentum. Here the orbitals with 
M = 0,1,2,3,4,5 are occupied, and the UHF solution is 
identical to the RHF solution with total angular momen- 
tum M tot = 15. 

Finally, in Fig. I14f e') for strong magnetic field we have 
a distinctly localized fivefold Wigner molecule. Figure 
I14f f) for Cj c — 2.5 shows a sixfold isomer which is higher 
in energy by 0.009 than the fivefold ground state. 



VI. CONCLUSION 

In conclusion, we have discussed the properties of un- 
restricted Hartree-Fock (UHF) calculations for electrons 
in a quantum dot, focusing on the regime of strong cor- 
relations, when the electrons begin to form a Wigner 
molecule. The UHF energies are good estimates of the 
true ground-state energies, especially for the polarized 
states, even at strong interaction. In this regime, the 
UHF energies become nearly spin independent, faster 
than it is the case for the true energies. However, the 
energy differences between different spin states cannot 
be resolved correctly by UHF, the polarized state is un- 
physically favored for stronger interaction. 

Regarding the interpretation of other quantities ob- 
tained from the UHF Slater determinant, we have shown 
that considerable caution must be taken: we find de- 
formed densities in the regime of intermediate interac- 
tion A ~ 1 ... 4. For stronger interaction the densities are 
azimuthally modulated for an even number of electrons 
per spatial shell, and round for an odd number per shell. 
The onset of this modulation is enhanced within UHF, so 
that UHF leads to an overestimation of the value of the 
critical density for the crossover to the Wigner molecule. 
We want to emphasize that the even-odd effect we found 
is an artifact of the symmetry breaking of UHF and arises 
from a degeneracy of states with opposite total angular 
momentum. 

For very strong interaction, we have shown that the 
UHF Hamiltonian corresponds to a tight-binding model 



of a particle hopping between the sites of the Wigner 
molecule. From the UHF orbital energies we have ob- 
tained the hopping matrix elements. This correspon- 
dence explains why the UHF energies become nearly spin 
independent which is expected for localized electrons and 
was not found with restricted HF. 

The maxima of the UHF densities mirror the classi- 
cal filling scheme with the electrons arranged in spatial 
shells. In contrast, the UHF two particle density (condi- 
tional probability density) has no direct physical mean- 
ing, because the UHF method can not take correlations 
properly into account. Finally, in a strong magnetic field 
the UHF densities are always molecule-like and there is 
no even-odd effect. 

The numerical complexity of the UHF method is com- 
parable to the frequently used density-functional ap- 
proach. However, as shown here, UHF has the advan- 
tage to cope also with the strongly interacting limit and 
gives further physical insight in that case. For the tiny 
energy differences which determine the spin ordering or 
the addition energies at A > 2 one has to employ the 
computationally more expensive quantum Monte Carlo 
methods. 
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APPENDIX A: CONFIGURATIONS OF 
CLASSICAL POINT CHARGES 

In TableE|we give the classical configurations for up to 
seven 2D electrons in a parabolic confinement potential 



N Geometry rl/X r s /A 4/3 E/r 2 a 

2 
3 
4 
5 
5* 
6 
6* 
7 



dumbbell (2) 


i 

4 




« 1.260 


3 


triangle (3) 


^« 0.577 


« 1.442 


9 
2 


square (4) 




: 0.957 


« 1.394 


6 


pentagon (5) 




« 1.376 


w 1.308 


15 
_> 


square (4,1) 




: 1.957 


« 1.251 


6 


pentagon (5,1) 


1 + y/l 




» 1.334 


15 

2 


hexagon (6) 


2 + J- a 


I 1.827 


« 1.223 


9 


hexagon (6,1) 


4 ^ v-3 


s 2.827 


« 1.414 


9 



12 




with zero magnetic field. r a is the distance of the outer 
electrons from the center measured in oscillator length Iq. 
r s is the nearest neighbor distance measured in effective 
Bohr radii ag. Energies are given in units of fko. These 
quantities depend only on N and A. 

For N — 5 and 6 we specify isomers with higher ener- 
gies. Due to the classical virial theorem there is a simple 
relationship between the energy and r a . When we denote 



the distance of the i-th electron from the center by r, , we 
have 

„ N 

E =\Y,rl- (Al) 

i=l 
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